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Abstract 

Two  cases  were  computationally  investigated  using  the  detaehed  eddy  simulation  (DES) 
turbulenee  model:  an  unmodified  Pak-B  blade  and  a  Pak-B  blade  with  a  dimple  loeated  at 
65%  of  axial  ehord.  Both  cases  were  created  so  that  they  simulated  an  infinite  span  with 
an  infinite  number  of  dimples.  The  eases  were  run  for  an  inlet  Reynolds  number  of 
25,000.  The  eomputed  results  were  used  to  resolve  the  loeation  of  separation  and 
reattaehment,  visualize  the  streamlines  for  the  dimpled  ease,  build  veloeity  magnitude 
eontour  and  veetor  plots,  and  map  the  thiekness  of  the  boundary  layer.  The  results  were 
then  eompared  to  previous  eomputational  and  experimental  studies  in  order  to  validate 
the  detached  eddy  simulation  model  for  future  researeh  into  the  effeet  of  dimples  on  low 
pressure  turbine  flow  fields.  Eor  the  unmodified  blade,  the  performanee  of  the  DES 
model  compared  favorably  to  other  available  viscous  and  turbulence  models.  Eor  the 
dimpled  blade,  preliminary  results  also  indieate  favorable  results  although  further 
development  of  the  flow  field  is  neeessary  to  verify  this.  Based  on  these  results,  future 
researehers  studying  dimples  on  turbine  blades  should  strongly  eonsider  using  the  DES 
turbulenee  model. 
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DETACHED  EDDY  SIMUEATION  ANAEYSIS  OF  PAK-B  EOW 
PRESSURE  TURBINE  BEADE 


CHAPTER  1.  INTRODUCTION 


Currently,  the  Department  of  Defense  uses  Unmanned  Aerial  Vehieles  (UAVs), 
like  the  Navy’s  X-47  Pegasus  pietured  in  Fig.  1,  for  several  missions  sueh  as  high- 
altitude  reeonnaissance  that  have  historically  been  performed  by  manned  aircraft.  As  the 
technology  advances,  UAVs  are  likely  to  take  on  an  ever-greater  number  of  these 
missions.  Unfortunately,  the  capability  of  UAVs  to  perform  these  missions,  which  often 
require  flight  at  high  altitudes  and  low  velocities,  is  severely  constrained  by  limitations 
inherent  in  low  pressure  turbines  (LPTs),  which  provide  power  to  the  aircraft. 
Specifically,  these  flight  conditions  can  result  in  an  operating  Reynolds  number  of  less 
than  25,000.  This  operating  condition  causes  the  development  of  a  laminar  boundary 
layer  over  the  turbine  blades.  The  laminar  nature  of  the  boundary  layer  combined  with 
the  large  turning  angles  generally  associated  with  LPTs  tends  to  cause  flow  separation 
over  much  of  the  turbine  blade  [Rouser,!].  This  phenomenon  limits  the  aircraft’s  range. 


1 


Figure  1;  X-47  Pegasus  UCAV  [2] 


altitude,  and  the  amount  of  eleetrieal  power  that  ean  be  extracted  from  the  engine.  Most 
current  research  emphasizes  work  on  delaying  or  eliminating  flow  separation  at  low 
Reynolds  numbers  without  adversely  affecting  blade  performance  at  high  Reynolds 
numbers. 

Engineers  at  Pratt  and  Whitney®  designed  the  Pak-B  blade  profile  specifically  to 
study  and  improve  the  performance  of  LPTs  at  low  Reynolds  numbers.  The  2-D  Pak-B 
shape  comes  from  the  3-D  Pak  blade  developed  for  commercial  use  at  Pratt  and  Whitney. 
In  order  to  analyze  the  Pak-B  shape  in  3-D  space,  the  basic  profile  is  extruded  in  the 
span-wise  direction.  However,  though  this  blade  design  is  useful  for  research  purposes,  it 
is  not  a  commercially  viable  design  since  it  neglects  taper  and  twist  [Rouser,!]. 

Lake  [3]  investigated  passive  flow  separation  control  techniques  on  a  Pak-B 
blade.  He  concluded  that  placing  dimples  on  the  blade  just  forward  of  the  natural  chord- 
wise  separation  location  tends  to  prevent  flow  separation.  Using  this  technique,  he  found 
that  losses  were  reduced  by  51.7%  for  a  blade  operating  at  a  Reynolds  number  of  45,000 
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and  freestream  turbulence  of  4%.  Lake  [3]  hypothesized  that  the  dimples  acted  as  vortex 
generators.  In  effect,  they  energize  the  flow  and  force  it  to  transition  from  laminar  to 
turbulent  before  it  reaches  its  natural  separation  point. 

There  are  many  advantages  to  passive  control  techniques  such  as  the  use  of 
dimples.  Unlike  other  methods,  they  have  been  shown  to  reduce  operating  inefficiencies 
at  low  Reynolds  numbers  and  they  do  not  adversely  affect  blade  performance  when 
operating  at  high  Reynolds  numbers  [Rouser,!;  Casey ,4].  Also,  they  neither  incur  a 
weight  penalty  when  compared  to  an  unmodified  blade  nor  do  they  require  energy 
beyond  the  basic  turbine  requirement  from  the  engine  in  order  to  operate.  Furthermore, 
they  can  be  retrofitted  into  existing  designs  without  extensive  redesign.  Due  to  these 
advantages,  the  use  of  dimples  for  separation  control  produces  favorable  results. 
However,  for  future  designs  to  effectively  make  use  of  this  technique,  a  more  complete 
understanding  of  the  flow  field  and  how  it  is  affected  by  dimples  is  required. 

The  current  research  builds  on  previous  studies  into  the  effect  of  dimples  on  low 
pressure  turbine  blades.  Two  blade  designs,  an  unmodified  Pak-B  blade  and  a  Pak-B 
blade  with  a  dimple  positioned  at  65%  of  axial  chord,  were  computationally  investigated 
using  a  recently  developed  turbulence  model  known  as  detached  eddy  simulation  (DBS). 

1.1  Low  Pressure  Turbine  Operation  at  Low  Reynolds  Numbers 

Unlike  other  turbomachinery  such  as  compressors  or  high  pressure  turbines,  low 
pressure  turbines  are  not  generally  subjected  to  extreme  environmental  conditions  that 
require  precautionary  measures,  such  as  active  cooling,  for  safe  operation.  Therefore, 
research  on  low  pressure  turbines  has  tended  to  focus  on  optimizing  their  aerodynamic 
characteristics  and  maximizing  their  loading  capabilities  and  efficiency.  Research  that 
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aims  to  improve  the  ability  of  low  pressure  turbines  to  operate  at  low  Reynolds  numbers 
falls  into  this  category. 

Low  pressure  turbines  always  operate  in  an  unsteady  environment.  As  they  rotate, 
the  turbine  blades  periodically  shed  wakes  that  generate  high  local  freestream  turbulence 
as  they  travel  downstream.  During  these  periods  of  high  turbulence,  the  boundary  layer 
generally  does  not  separate  from  the  blade.  However,  once  a  wake  has  passed  and  the 
turbulence  decreases,  the  flow  again  becomes  laminar  and,  as  a  result,  more  vulnerable  to 
separation  [Lake,2]. 

1.2  Previous  Experiments  in  Flow  Separation  on  LPT  Blades 

Numerous  experiments  dealing  with  low  pressure  turbines  operating  at  Reynolds 

numbers  below  200,000  have  greatly  enhanced  our  knowledge  of  the  aerodynamics  at 
these  conditions.  Sharma  et  al.  [5]  found  that  the  loss  coefficient  increased  dramatically 
when  the  operating  Reynolds  number  dropped  below  95,000.  Murawski  et  al.  [6]  and 
Qiu  and  Simon  [7]  subsequently  showed  that  operating  at  low  Reynolds  numbers 
significantly  decreased  the  effectiveness  of  low  pressure  turbine  blades.  Specifically, 
Murawski  et  al.  [6]  noted  separation  from  75%  to  90%  of  axial  chord  at  Reynolds 
numbers  as  low  as  53,000  and  Qiu  and  Simon  [7]  verified  that  reattachment  occurs 
downstream  of  90%  of  axial  chord.  Simon  and  Volino  [8]  ran  experiments  on  a  Pak-B 
blade  operating  at  flow  conditions  similar  to  those  experienced  by  high-altitude 
reconnaissance  aircraft  and  confirmed  the  decrease  in  effectiveness. 

These  research  efforts  created  a  foundation  for  Lake’s  [3]  research  with  a  Pak-B 
blade  operating  at  low  Reynolds  numbers.  Lake  [3]  showed  that  several  indicators  of 
effectiveness,  including  the  total  pressure  loss  coefficient  and  wake  momentum  deficits. 


4 


increased  due  to  flow  separation  at  these  operating  eonditions.  At  a  Reynolds  number  of 
45,000,  he  discovered  a  natural  separation  loeation  at  70%  of  axial  ehord  [Lake,3]. 
Subsequently,  he  attempted  to  control  separation  with  various  deviees,  ineluding  trip 
wires,  span-wise  V-grooves,  and  reeessed  dimples.  Though  he  noted  a  reduetion  in  flow 
separation  with  eaeh  of  these  devices,  he  eoncluded  that  dimples  most  effeetively  redueed 
separation. 

Rouser  [1]  further  deereased  the  operating  Reynolds  number  to  25,000  and 
experimentally  investigated  several  other  dimples  patterns  eonsisting  of  either  full  or 
stream-wise  half-dimples  loeated  at  50%,  55%,  or  65%  of  axial  ehord.  Of  these 
eonfigurations,  he  found  the  most  effeetive  to  be  a  row  of  full  dimples  loeated  at  65%  of 
axial  chord.  This  design  redueed  the  amount  of  separation  by  28%.  He  also  studied 
eomputational  models  of  both  an  unmodified  Pak-B  blade  and  a  Pak-B  blade  modified 
with  both  full  and  half-dimples  loeated  at  65%  of  axial  ehord  in  order  to  gain  more 
knowledge  of  the  flow  structure.  In  his  eomputational  research,  he  utilized  several 
viseous  models,  ineluding  a  laminar  solver  and  the  Spallart-Allmaras,  K-m,  and  K-s 
turbulence  models.  Rouser’s  [1]  eomputational  results  described  a  highly  transitional 
flow  field.  Unfortunately,  the  turbulenee  models  amplified  the  instability  inherent  to  the 
flow  field  so  that  they  eould  not  resolve  a  separation  loeation.  The  laminar  solver 
resolved  the  separation  loeation  at  66%  of  axial  ehord. 

Casey  [4]  expanded  the  base  of  research  by  investigating  several  new  dimple 
patterns.  Experimentally,  Casey  [4]  studied  a  row  of  dimples  with  a  2.22  em  eenter-to- 
eenter  spacing  located  at  65%  of  axial  chord,  a  row  of  dimples  with  a  4.44  cm  center-to- 
center  spacing  located  at  65%  of  axial  ehord,  and  a  two-row  staggered  pattern  with 
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dimples  located  at  65%  and  76%  of  axial  chord  with  a  4.44  cm  center-to-center  spacing. 
The  results  from  each  dimple  pattern  were  compared  to  results  from  an  unmodified  blade 
at  operating  Reynolds  numbers  of  25,000,  45,000,  and  100,000.  He  found  that  no  dimple 
pattern  performed  significantly  more  effectively  than  any  other  dimple  pattern  and  that 
each  pattern  reduced  the  average  total  pressure  loss  coefficient  by  34%  for  an  operating 
Reynolds  number  of  25,000  and  a  freestream  turbulence  of  1%.  Casey  [4]  also  created 
slightly  modified  computational  models  of  these  cases  and  solved  each  case  using  a 
laminar  solver.  These  models  allowed  Casey  [4]  to  create  a  visualization  of  the  flow  field. 
Comparisons  to  the  experimental  data  verified  the  accuracy  of  the  computational  results. 

1.3  Current  Research  Objectives 

The  objective  of  the  current  research  is  to  investigate  the  ability  of  a  numerical 
model  called  detached  eddy  simulation  (DBS)  to  resolve  the  flow  field  surrounding  a  low 
pressure  turbine  blade.  In  order  to  allow  valid  comparisons  with  Casey’s  [4]  results,  two 
blade  geometries  were  analyzed:  an  unmodified  Pak-B  blade  and  a  Pak-B  blade  with  a 
spherical  dimple  located  at  65%  of  axial  chord.  A  unique  grid  was  created  for  both  of 
these  blade  geometries.  Both  grids  employed  periodic  boundary  conditions  in  both  the 
span-wise  and  pitch-wise  directions  in  order  to  simulate  an  infinitely  long  2-D  blade  with 
a  full  span  of  dimples  and  an  infinite  linear  cascade,  respectively.  Both  grids  were  built 
in  Gridgen®  and  solved  in  Fluent®  on  a  time-accurate  basis  using  the  DBS  turbulence 
model.  Post  processing,  including  calculation  of  separation  and  reattachment  locations, 
streamline  visualization,  velocity  magnitude  contour  plots,  and  boundary  layer  profiles, 
was  performed  in  Fieldview®. 
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1.4  Chapter  Summary 

Chapter  1  contains  an  introduction  to  the  current  research  and  some  background 
regarding  its  origins.  Chapter  2  discusses  the  theory  behind  the  analysis  of  low  pressure 
turbines,  various  methods  used  in  attempts  to  suppress  flow  separation,  and  the  theory 
and  application  of  the  DBS  turbulence  model.  Chapter  3  discusses  the  computational 
methodology  and  results.  Chapter  4  sums  up  the  conclusions  of  this  thesis  and  provides 
recommendations  for  future  work  on  this  topic. 
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CHAPTER  2.  BACKGROUND  AND  THEORY 


The  primary  focus  of  this  research  is  the  advancement  of  the  technology  found  in 
UAVs.  Specifically,  it  is  desired  to  enhance  the  operational  capacities  of  these  aircraft  by 
optimizing  parameters  such  as  range,  altitude,  and  endurance.  In  order  to  do  so,  methods 
by  which  turbine  efficiencies  can  be  increased  were  investigated.  Pressure  and  velocity 
are  the  computational  parameters  most  directly  correlated  with  turbine  efficiency. 
Therefore,  pressure  and  velocity  data  were  gathered  from  the  simulations  and  used  during 
data  analysis. 

The  operational  conditions  experienced  by  the  low  pressure  turbines  used  to 
operate  UAVs  often  reach  Reynolds  numbers  as  low  as  25,000.  Previous  research  and 
experience  has  shown  that  the  flow  over  the  turbine  blade  separates  under  these 
conditions.  Since  it  is  believed  that  the  efficiency  of  the  turbine  will  increase  if  the 
amount  of  flow  separation  can  be  reduced  or  eliminated,  this  research  investigates 
methods  to  minimize  flow  separation.  These  research  goals  require  a  detailed  model  of 
the  behavior  of  the  boundary  layer.  The  computational  investigation  during  the  course  of 
this  research  provides  a  model  of  the  boundary  layer  as  well  as  a  picture  of  the  rest  of  the 
flow  field. 


2.1  Low  Pressure  Turbine  Performance  Evaluation 

The  maximum  range  of  an  aircraft  is  often  used  as  a  measure  of  its  engine 
performance.  Breguet’s  range  formula  shows  the  inversely  proportional  relationship 
between  aircraft  range  and  the  engine’s  thrust  specific  fuel  consumption  (TSFC): 


C 

R  =  -  ^  In 
Cr 


W 


final 


\  '^initial  J 


U  1 


TSFC  g 


(1) 
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where  U  is  the  flight  veloeity,  g  is  the  aeceleration  due  to  gravity,  Cd  is  the  drag 
coefficient,  Cl  is  the  lift  coefficient,  R  is  the  aircraft  range,  Wfmai  is  the  aircraft’s  final 
weight,  and  Wmitiai  is  the  aircraft’s  initial  weight  [Mattingly,9;46].  The  TSFC  is  given  by: 

TSFC  =  (2) 

nohpr 

where  rjo  is  the  engine  efficiency  and  hpr  is  the  lower  heating  value  of  the  fuel  used  by  the 
engine  [Mattingly,  11:30].  Substituting  this  relationship  into  Breguet’s  range  formula 
gives: 


R  = 


(3) 


These  algebraic  manipulations  show  that  as  the  efficiency  of  the  engine  increases,  the 
maximum  range  of  the  aircraft  also  increases.  Since  the  increase  in  engine  efficiency  also 
reduces  the  amount  of  fuel  required  for  any  arbitrary  range,  improvements  in  engine 
efficiency  can  simultaneously  increase  aircraft  mission  effectiveness  and  reduce  the 
amount  of  money  that  would  otherwise  need  to  be  spent  on  fuel  for  the  aircraft. 

The  Brayton  cycle  describes  an  ideal  thermodynamic  cycle  for  a  gas  turbine 
engine.  See  Fig.  2  for  a  graphical  representation  of  the  Brayton  cycle  and  its  associated 
T-S  diagram,  where  T  is  the  temperature  and  S  is  the  entropy. 
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Figure  2;  T-S  Diagram  and  Illustration  of  Brayton  Cycle  [Mattingly,  9:233] 


Between  stations  0  and  3,  the  working  fluid  undergoes  isentropic  compression,  resulting 
in  an  increase  in  both  temperature  and  pressure,  with  the  compressor  mechanism  located 
between  stations  2  and  3.  Between  stations  3  and  4,  it  is  then  heated  by  the  combustion 
of  fuel  injected  into  a  combustor.  This  process  takes  place  at  a  constant  pressure  and  is 
often  more  simply  modeled  as  a  heat  input  as  shown  in  Fig.  2.  Then,  between  stations  4 
and  9,  the  working  fluid  passes  through  a  turbine  and  a  nozzle,  resulting  in  an  isentropic 
decrease  in  both  temperature  and  pressure.  Then,  between  stations  9  and  0,  the  working 
fluid  passes  through  another  heat  exchanger  and  rejects  heat  from  the  system.  After 
passing  through  every  station,  the  working  fluid  now  has  the  same  temperature,  pressure, 
and  entropy  as  it  had  before  the  cycle  started  [Oates,  10:136]. 

The  first  law  of  thermodynamics  states  that  for  a  closed  system: 

^E  =  Q-W  (4) 

where  AC  is  the  net  change  in  energy,  Q  is  the  heat  added  to  or  removed  from  the  system, 
and  W  is  the  net  work  transferred  into  or  out  of  the  system  [Moran  and  Shapiro,  1 1:58]. 
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Assuming  the  working  fluid  is  a  calorically  perfect  gas,  analysis  of  the  Brayton  cycle 
yields  the  following  relationships: 


Wc=mc^(T,,-TJ  (5) 

W,=mc^(T,,-TJ  (6) 


where  W  is  the  rate  that  work  passes  through  the  system,  m  is  the  mass  flow  rate,  Cp  is  the 
average  specific  heat  at  constant  pressure,  and  Tt  is  the  stagnation  temperature 
[Mattingly,  1 1 :234].  In  the  Brayton  cycle,  the  work  required  by  the  compressor  is 


supplied  by  the  turbine  so  that  =W^  . 

As  stated  previously,  the  Brayton  cycle  represents  an  ideal  thermodynamic  cycle. 
As  such,  it  ignores  the  effects  of  entropy  production.  All  real  systems  must  also  follow 
the  second  law  of  thermodynamics,  which  states  that  for  a  closed  system: 


ds  > 


dQ 

T 


(V) 


where  ds  is  the  net  change  in  entropy  and  dQ  is  the  net  change  in  the  amount  of  heat  in 
the  system  [Anderson,  12:27].  Put  simply,  any  real  process  will  cause  the  entropy  of  the 
system  and  its  surroundings  to  increase.  Therefore,  the  Brayton  cycle  can  only  give  an 
approximate  representation  of  a  gas  turbine  engine.  A  real  turbine  engine  experiences 
losses  due  to  many  factors  including  flow  separation,  a  potentially  major  factor  in  the 
overall  engine  efficiency. 

A  T-S  diagram  can  illustrate  a  non-ideal  system  in  the  same  manner  that  it 
described  the  Brayton  cycle  above.  Fig.  3  shows  a  T-S  diagram  for  a  non-ideal  turbine 
stage. 
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Figure  3;  T-S  Diagram  for  a  Non-ideal  Turbine  Stage  [Lake,  3] 


As  seen  in  Fig.  3,  entropy  inereases  over  both  the  stator  and  the  rotor.  Furthermore,  the 
pressure  decreases  over  the  stator  and  increases  over  the  rotor.  The  Gibbs  function 
relates  the  change  in  entropy  to  the  changes  in  temperature  and  pressure: 


ds  =  Cp 


(8) 


where  ds  is  the  change  in  entropy,  R  is  the  gas  constant,  Cp  is  the  average  specific  heat  at 
constant  pressure,  Tt  is  the  stagnation  temperature,  dTt  is  the  change  in  stagnation 
temperature,  Pt  is  the  stagnation  pressure,  and  dPt  is  the  change  in  stagnation  pressure 
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[Oates,  10:38].  As  seen  in  Fig.  3,  the  stagnation  temperature  is  constant  across  the  stator. 
Neglecting  this  term  in  Eq.  8  gives: 


ds  =  -R 


P, 


(9) 


Eq.  9  describes  the  relationship  between  the  pressure  drop  and  the  amount  of  entropy 
production  across  the  stator.  Eogically,  since  there  are  relationships  between  pressure 
drop  and  entropy  production  and  between  entropy  production  and  stage  efficiency,  it  is 
also  possible  to  directly  relate  pressure  drop  and  stage  efficiency: 


V, 


1-E 


1 

\-7T,  /r. 


(10) 


where  r|t  is  the  turbine  efficiency,  /t  is  the  ratio  of  specific  heats  of  the  gases  within  the 
T  P 

turbine,  r,  =  — ,  and  n,  =  —  [Oates,  221].  Assuming  a  fixed  stagnation  temperature 

Pt\  Pt\ 


ratio,  decreasing  the  exit  stagnation  pressure,  which  can  be  measured  within  the  linear 
cascade  of  a  gas  turbine  engine,  also  decreases  the  turbine  efficiency. 

The  total  pressure  loss  coefficient,  another  measure  of  efficiency,  is  defined  by 
the  following  equation: 


CO  = 


P  -P 

t, inlet  t.exit 


^  inlet 


(11) 


where  ro  is  the  total  pressure  loss  coefficient,  Pt  is  the  stagnation  pressure,  and  qmiet  is  the 
dynamic  pressure  at  the  inlet  [Casey, 4].  As  shown  by  Eq.  1 1,  an  increase  in  the  exit 
stagnation  pressure  will  decrease  the  total  pressure  loss  coefficient.  Since  minimizing  or 
eliminating  the  flow  separation  will  increase  the  exit  stagnation  pressure,  it  will  also 
decrease  the  total  pressure  loss  coefficient.  Eake  [3]  showed  this  effect  in  Eig.  4. 
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Figure  4;  T-S  Diagram  Showing  Improvements  Due  to  Higher  Exit  Pressures  [Lake, 3] 

The  baseline  plot  represents  a  turbine  with  no  surface  modifications  and  the  improved 
plot  shows  the  significant  improvement  in  efficiency  that  can  result  from  only  a  small 
increase  in  the  exit  stagnation  pressure. 

In  his  research,  Lake  [3]  found  that  dimples  etched  into  the  surface  of  a  Pak-B 
blade  reduced  the  total  pressure  loss  coefficient.  Later,  Rouser  [1]  confirmed  this  result 
and  showed  that  a  dimple  located  at  65%  of  the  axial  chord  had  the  greatest  effect  on 
increasing  the  engine  efficiency.  Casey  [4]  investigated  several  dimple  patterns  and 
found  that  each  performed  to  a  similar  degree.  Ultimately,  he  concluded  that  a  single  row 
of  dimples  at  65%  of  axial  chord  spaced  4.44  cm  from  each  other  produced  the  best 
results  at  the  lowest  cost  [Casey, 4]. 
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2.2  Low  Pressure  Turbine  Boundary  Layer  Aerodynamics 

The  flow  fields  about  low  pressure  turbine  blades  are  inherently  unsteady  and 
transitional.  Halstead  et  al.  [13;  14]  showed  that  turbulence  levels  between  1%  and  20% 
exist  in  the  freestream.  They  also  detailed  other  structures  within  the  flow  field  such  as 
passing  wakes,  centrifugal  effects,  and  vortices.  Various  researchers  have  shown  that  the 
boundary  layers  on  low  pressure  turbine  blades  operating  at  low  Reynolds  numbers 
generally  transition  from  laminar  to  turbulent  over  the  suction  surface  of  the  blade 
[Mayle,15;  Addison  and  Hodson,  16].  This  transition  is  almost  steady  but  is  affected  by 
time-dependent  wakes  caused  by  the  relative  motion  of  the  rotor  and  the  stator.  As  these 
wakes  move  downstream,  the  regions  through  which  they  pass  quickly  switch  from 
transitional  to  turbulent  and  back  to  transitional  once  the  wakes  have  moved  on. 

2.2,1  Experimental  Laminar-to-Turbulent  Boundary  Layer  Transition 

Many  researchers  have  focused  their  work  on  efforts  to  determine  the  exact 
location  of  transition  and  the  length  of  the  blade  surface  under  the  transitional  region  in 
order  to  determine  a  relationship  between  the  duration  of  transition  and  performance 
losses.  These  variables  are  determined  primarily  by  the  transitional  modes  of  the  flow 
field.  Walker  [17]  identified  these  modes  as  natural  transition,  bypass  transition,  and 
separation  bubble  transition.  Unfortunately,  since  many  factors  can  affect  these 
transitional  modes,  including  Reynolds  number,  freestream  turbulence,  and  pressure 
gradients,  research  on  this  topic  has  been  slow  and  difficult. 

Each  of  these  modes  exists  under  very  different  conditions.  Natural  transition 
only  occurs  in  flow  fields  characterized  by  low  freestream  turbulence  and  the  absence  of 
adverse  pressure  gradients.  This  process  is  initialized  by  the  formation  of  two- 
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dimensional  Tollmien-Schlichting  (T-S)  waves.  As  the  flow  progresses  and  the  T-S 
waves  beeome  unstable,  they  ereate  three-dimensional  loop  vortiees  and  large 
fluetuations  [Mayle,  15].  These  struetures  merge  and  beeome  turbulent  spots  that 
inerease  in  size  as  they  move  downstream  until  they  fill  the  entire  boundary  layer.  At 
that  point,  the  flow  is  eonsidered  fully  turbulent. 

Flow  fields  with  freestream  turbulenee  levels  of  approximately  20%  or  higher 
tend  to  experienee  bypass  transition.  T-S  waves  do  not  form  in  these  flow  fields  but  the 
higher  level  of  turbulenee  still  allows  the  formation  of  turbulent  spots  similar  to  those 
formed  by  natural  transition.  Walker  [17]  showed  that  linear  stability  theory  eould  still 
be  used  to  model  these  flow  fields  and  the  disturbanees  within  them  so  that  the  length  of 
the  transitional  region  ean  be  determined. 

A  laminar  flow  field  eharaeterized  by  adverse  pressure  gradients  often  sees  the 
development  of  separation  bubbles  in  whieh  transition  oeeurs  in  the  shear  layer  between 
the  freestream  and  the  separation  bubble.  Separate  investigations  by  Mayle  [14]  and  by 
Qiu  and  Simon  [8]  eonfirmed  that  reattaehment  of  the  boundary  layer  to  the  blade  surfaee 
is  eaused  by  the  transition  of  the  flow  field  from  laminar  to  fully  turbulent.  This 
transition  mode  eharaeterizes  the  flow  field  of  a  Pak-B  blade  operating  at  low  Reynolds 
numbers. 

2,2.2  Numerical  Models  for  Low  Pressure  Turbine  Performance 

Previous  numerieal  studies  of  low  pressure  turbines  have  been  hindered  by  the 
inability  of  the  models  used  to  resolve  the  flow  field.  This  diffieulty  has  been  attributed 
to  the  absenee  of  experimental  researeh  to  eompare  with  those  numerieal  models  [Werle, 
18].  Halstead  et  al.  [19]  attempted  to  use  many  different  models,  ineluding  STANX, 
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KEP,  Navier-Stokes,  and  Fan-Lakshminarayana,  that  ultimately  failed  to  produee  realistie 
results.  Lake  [3]  used  the  two-dimensional  Allison  Vane -Blade  Interaetion  (VBI)  eode. 
He  found  that  though  this  code  generally  worked  well,  it  still  had  difficulties  resolving 
the  separated  region  at  low  Reynolds  numbers  and  was  unable  to  deal  with  three- 
dimensional  flow  structures  characteristic  of  the  flow  fields  surrounding  dimples.  As 
such,  he  used  this  code  for  an  unmodified  Pak-B  blade  only. 

For  his  research.  Lake’s  [3]  computational  model  included  an  0-grid  around  the 
blade  and  an  H-grid  in  the  blade’s  passage.  Rouser  [1]  built  on  Lake’s  [3]  model  in  order 
to  create  a  three-dimensional  grid.  He  extruded  a  structured  0-grid  from  the  surface  of 
the  blade  to  produce  hexahedral  cells  in  the  near-wall  region.  He  then  modeled  the  blade 
passage  using  an  unstructured  grid  characterized  by  tetrahedral  cells.  The  flow  region 
directly  above  the  dimples  was  also  modeled  using  an  unstructured  grid  style.  However, 
this  region  used  prism  cells  instead  of  hexahedral  cells  and  extended  all  the  way  to  the 
wall.  See  Fig.  7  for  an  illustration  of  these  cell  types. 
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Figure  5;  Illustration  of  Cell  Types  [Casey,4] 
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After  building  the  grid,  Rouser  [1]  computed  the  flow  field  in  Fluent®.  Initially, 
he  assumed  a  turbulent  flow  and  employed  a  Reynolds  Averaged  Navier-Stokes  (RANS) 
turbulence  model  throughout  the  domain.  He  found  that  the  RANS  model  predicted  that 
no  separation  would  occur  in  a  flow  field  with  an  inlet  turbulence  level  of  1%.  After 
developing  this  solution,  Rouser  [1]  switched  to  a  laminar  solver  and  discovered  that  it 
could  properly  resolve  the  flow  field  at  the  separation  and  reattachment  locations  but  was 
unable  to  do  so  in  the  turbulent  region  aft  of  the  reattachment  point.  He  surmised  that  a 
new  turbulence  model  able  to  resolve  this  region  properly  was  necessary. 

Casey  [4]  used  the  same  overall  grid  style  as  Rouser  [1].  However,  rather  than 
attempt  to  use  a  turbulence  model,  he  computed  the  flow  field  using  a  time-accurate 
laminar  solver  since  he  believed  this  to  be  more  physically  accurate  than  running  the 
laminar  solver  from  a  steady  turbulent  solution  as  Rouser  [1]  had  done.  This  approach 
was  more  successful  than  previous  attempts  to  describe  the  flow  field.  Unlike  Rouser  [1], 
Casey  [4]  resolved  large  vortical  structures  in  the  flow  field  with  a  dimpled  blade. 

Rizzeta  and  Visbal  [20]  performed  a  computational  investigation  using  the  large- 
eddy  simulation  (LES)  turbulence  model  on  an  unmodified  Pak-B  blade  with  operating 
Reynolds  numbers  as  low  as  25,000.  They  found  that  the  LES  turbulence  model  was  able 
to  capture  the  features  of  the  transitional  flow  field.  Also,  their  results,  which  included 
variables  such  as  coefficient  of  pressure,  mean  velocity,  and  fluctuating  velocity, 
compared  fairly  well  with  a  two-dimensional  model  of  the  blade.  However,  the  LES 
model  did  not  provide  perfect  results.  The  coefficient  of  pressure  data  differed  from 
experimental  data  although  they  attributed  this  error  to  the  non-periodic  nature  of  the 
experimental  data.  Prasad  et  al.  [21]  used  the  LES  model  to  resolve  the  flow  field  for  an 
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HEV  static  mixer  and  discovered  that,  whereas  previous  turbulence  models  could  not 
resolve  experimentally  observed  vortices,  the  LES  model  predicted  their  existence. 

The  current  research  employed  a  turbulence  model  called  detached  eddy 
simulation  (DES)  throughout  the  flow  field.  DES  is  a  hybrid  model  that  applies  different 
numerical  turbulence  models  to  different  regions  of  the  grid.  Specifically,  it  uses  a 
RANS  turbulence  model  in  the  boundary  layer  and  a  large-eddy  simulation  approach 
outside  of  the  boundary  layer.  The  DES  model  used  in  this  research  is  based  on  the 
Spallart-Allmaras  (S-A)  turbulence  model.  In  order  to  modify  the  S-A  model,  the  DES 
model  replaces  the  S-A  variable  dw,  the  distance  to  the  nearest  wall,  with: 

d  =  rmn{d^,CDEs^)  (12) 

where  Cdes  is  the  DES  model  constant  and  A  is  the  largest  distance  between  the  cell 
center  being  considered  and  the  cell  center  of  its  neighbors  [Kapadia  et  ah,  22].  Inside 

the  boundary  layer,  r/  =  dw  and  the  normal  S-A  model  operates  so  that  RANS  calculations 
are  made.  Otherwise,  outside  the  boundary  layer,  the  model  becomes  grid-scale 
dependent  and  switches  to  the  EES  turbulence  model  [Kapadia  et  ah,  22]. 

Due  to  the  relative  newness  of  the  DES  model,  gridding  practices  are  still 
somewhat  underdeveloped.  Spalart  [23]  provides  fairly  complete  instructions  and  defines 
three  specific  flow  regions  within  any  grid  that  need  to  be  gridded  appropriately  in  order 
to  obtain  the  most  physically  accurate  results  in  the  most  cost-effective  manner:  the  Euler 
region,  the  RANS  region,  and  the  EES  region.  The  Euler  region  is  located  upstream  and 
to  the  sides  of  the  object  being  modeled  and  contains  no  turbulent  flow.  This  region 
should  generally  have  the  greatest  volume  but  the  least  number  of  grid  points  and  is 
marked  by  even,  isotropic  grid  spacing.  The  RANS  region  contains  the  boundary  layer 
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including  the  initial  location  of  separation  and  consists  of  inner  viscous  sub-region  and  an 
outer  sub-region.  The  inner  sub-region  should  be  modeled  in  the  same  manner  as  a  full- 
domain  RANS.  In  the  outer  sub-region,  the  grid  spacing  is  allowed  to  expand  such  that 
the  spacing  ultimately  equals  that  of  the  adjoining  Euler  region.  As  such,  this  region 
fulfills  the  needs  of  numerical  robustness  rather  than  any  physical  requirement.  The  LES 
region  is  the  primary  region  containing  turbulence  and  consists  of  three  sub-regions:  the 
viscous  sub-region,  the  focus  sub-region,  and  the  departure  sub-region.  The  viscous  sub- 
region  is  modeled  in  the  same  manner  as  the  viscous  sub-region  of  the  RANS  region. 

The  focus  region  sits  close  to  the  object  being  modeled  and  it  is  important  that  the 
turbulence  be  accurately  resolved  here  so  fairly  tight  grid  spacing  is  required.  In  the 
departure  region,  the  spacing  expands  significantly  since  the  structure  of  the  turbulence  is 
no  longer  as  important. 

2.3  Separation  Control  on  Low  Pressure  Turbines 

Research  by  Murawski  et  al.  [6]  and  Qiu  and  Simon  [7]  has  demonstrated  that 
Reynolds  numbers  below  300,000  or  freestream  turbulence  below  10%  more  easily 
allows  separation  than  operating  conditions  above  these  values.  Generally,  both 
conditions  are  met  by  the  mission  profile  of  a  high-altitude  reconnaissance  aircraft.  As 
the  Reynolds  number  decreases,  the  transition  location  on  the  blade  moves  downstream 
until  it  is  far  enough  that  the  flow  field  cannot  follow  the  curvature  of  the  blade  and 
separation  occurs.  Many  researchers  have  investigated  methods  by  which  low  pressure 
turbine  blades  can  be  modified  to  minimize  the  amount  of  separation. 
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2,3.1  Active  and  Passive  Flow  Separation  Control  Methods 


Separation  control  methods  can  be  classified  as  either  active  or  passive.  The 
operation  of  active  controls  requires  an  input  of  energy  whereas  passive  controls  do  not. 
Rather,  passive  control  methods  generally  involve  modifications  to  the  surface  of  the 
turbine  blade.  On  the  other  hand,  active  controls  can  be  switched  on  or  off  or  adjusted 
for  intensity  as  conditions  require  whereas,  once  instituted,  passive  controls  cannot  be 
adjusted.  Ultimately,  the  advantages  of  active  controls  make  them  better  suited  to  flow 
fields  with  varying  conditions  while  passive  controls  are  better  suited  to  flow  fields  with 
fairly  steady  operating  conditions. 

Johnston  and  Nishi  [24]  and  Compton  and  Johnston  [25]  showed  that  vortex 
generator  jets  (VGJs)  control  separation  on  fiat  plates  subjected  to  adverse  pressure 
gradients.  Bons  at  al.  [26]  investigated  the  interactions  between  VGJs  and  a  Pak-B  blade 
and  found  that  pulsating  VGJs  operate  as  effectively  as  steady  VGJs  with  less  mass  flow. 
Unfortunately,  this  type  of  active  control  requires  the  costly  input  of  high  pressure  air 
between  the  turbine  and  the  compressor.  Furthermore,  they  add  weight  to  the  engine  and 
can  easily  become  clogged  by  combustor  exhaust  particles.  As  such,  they  require  more 
maintenance  to  remain  effective.  Borgeson  [27]  attempted  to  use  micro¬ 
electromechanical  systems  (MEMS)  to  delay  separation  over  a  cylinder.  Though  he  met 
with  some  success,  this  method  cannot  currently  be  effectively  used  with  a  low  pressure 
turbine  due  to  the  inability  of  MEMS  to  withstand  extreme  operating  conditions  and  the 
difficulty  associated  with  installing  these  devices  in  a  low  pressure  turbine  blade. 

Passive  controls  tend  to  be  cheaper,  easier  to  use,  and  better  able  to  withstand 
harsh  operating  conditions.  Eurthermore,  by  improving  the  ability  of  current  blade 


21 


designs  to  avoid  separation,  they  prevent  eostly  redesigns  that  ean  adversely  affeet  the 
effleieney  of  low  pressure  turbine  blades  in  other  ways. 

The  most  eommon  passive  eontrol  teehniques  are  modifieations  to  the  surfaee  of 
the  low  pressure  turbine  blade.  Lin  et  al.  [28]  found  that  large-eddy  breakup  (LEBU) 
deviees  redueed  separation  at  small  positive  angles  of  attaek.  They  also  showed  that 
vortex  generators  submerged  within  the  blade  effeetively  suppressed  separation  so  long 
as  the  deviee  had  the  eorreet  height  and  was  properly  positioned  relative  to  the  natural 
separation  loeation  of  the  blade  surfaee.  Austin  [29]  and  Gamerdinger  and  Shreeve  [30] 
researehed  triangular  plow  vortex  generators  and  found  that  they  redueed  shoek-indueed 
separation.  Unfortunately,  these  deviees  also  redueed  downstream  stagnation  pressures. 
Rao  and  Kariya  [31]  looked  at  submerged  semi-eireular  and  simple  vane  vortex 
generators.  They  found  that  the  submerged  deviees  performed  better  than  the  vane 
deviees  as  long  as  their  design  ineluded  the  eorreet  height,  sweep  angle,  spaeing,  and 
length.  See  Fig.  6  for  illustrations  of  these  passive  eontrol  deviees. 


a)  Plow  h)  Semi-Ciaiilar  c)SnnplcVanc 

Figure  6:  Surfaee  Vortex  Generator  Deviees  [Fake,3] 
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Figure  7:  More  Vortex  Generator  Deviees  [Lake, 3] 


Due  to  the  interactions  between  the  flow  field  and  the  protruding  structures 
discussed,  these  devices  incur  a  drag  penalty  that  offsets  many  of  the  performance 
enhancements  gained  by  the  reduced  separation.  Lin  et  al.  [32]  proposed  several 
alternatives  to  these  protruding  structures  designed  to  reduce  the  drag  penalty  including 
a±  45°  small-scale  riblet  pattern  and  transverse  grooves.  The  transverse  grooves 
successfully  reduced  separation  while  the  riblet  pattern  actually  increased  the  amount  of 
separation.  McCormick  [33]  proposed  a  passive  cavity  method  that  uses  the  adverse 
pressure  gradient  behind  the  separation  location  to  assist  in  reducing  separation.  He 
found  that  this  idea  reduced  the  stagnation  pressure  loss  but  contributed  to  the 
deterioration  of  the  boundary  layer. 
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Other  separation  eontrol  ideas  inelude  leading  edge  vane  vortex  generators,  small- 
seale  V-groove  riblets,  small-seale  square  ribs,  and  a  vortex  generator  shaped  like  a 
NACA  64-812  airfoil  [Casey, 4],  Unfortunately,  all  of  these  teehniques  are  very  diffieult 
to  use.  Bloeh  and  Mueller  [34]  investigated  roughness  strips  plaeed  near  the  leading  edge 
of  an  airfoil  that  are  mueh  easier  to  install.  Lake  [3]  examined  three  passive  eontrol 
teehniques:  a  V-groove  ehannel,  a  boundary  layer  trip  wire,  and  dimples.  He  found 
dimples  to  be  the  most  effeetive  separation  suppression  teehnique  but  also  found  that 
their  effeetiveness  was  highly  dependent  on  several  faetors  ineluding  the  position  of  the 
dimple  relative  to  the  separation  loeation.  Rouser  [1]  and  Casey  [4]  both  eonfirmed  this 
eonelusion  and  eoneluded  that  the  best  loeation  for  the  dimple  was  at  65%  of  the  axial 
ehord. 
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Figure  8:  More  Vortex  Generator  Deviees  [Lake,3] 
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2,3.2  Previous  Research  in  Dimples  for  Boundary  Layer  Control 

Both  submerged  and  protruding  vortex  generators  have  been  shown  to  suppress 
boundary  layer  separation.  Protruding  surface  modifications  entrain  air  from  the 
freestream  into  the  boundary  layer,  thereby  energizing  it  and  reducing  the  amount  of 
separation  [Taylor,  35].  Also,  vortices  shed  off  the  leading  edge  which  mixes  the 
freestream  and  the  boundary  layer  downstream.  On  the  other  hand,  submerged  vortex 
generators  do  not  incur  as  severe  a  drag  penalty  as  the  protruding  modifications. 
Advantageously,  they  are  easier  to  install  on  low  pressure  turbine  blades  and  can  handle 
harsher  environmental  conditions.  But,  they  are  less  effective  unless  placed  very  near  the 
separation  location.  Lake  [3],  Rouser  [1],  and  Casey  [4]  have  all  shown  the  strong 
potential  that  dimples  have  in  boundary  layer  control  applications. 

Research  has  generally  confirmed  that  recessed  dimples  reduce  the  amount  of 
separation  in  the  boundary  layer.  Therefore,  more  recent  research  has  focused  more  on 
the  factors  that  can  affect  the  effectiveness  of  dimples  more  than  the  viability  of  dimples 
themselves.  Bearman  and  Harvey  [38]  investigated  spherical  and  hexagonal  dimples  on  a 
sphere  and  demonstrated  that  both  shapes  generated  downstream  vortices.  Of  the  two, 
they  found  that  the  hexagonal  dimple  shape  operated  more  efficiently. 

In  his  research.  Lake  [3]  built  on  the  procedures  and  results  of  Bearman  and 
Harvey  [38].  He  milled  dimples  into  a  Pak-B  blade  with  the  same  roughness 
characteristics  as  Bearman  and  Harvey  [38]had  previously  used  (k/D  =  9  x  10'  where  k 
is  the  maximum  dimple  depth  and  D  is  the  spherical  diameter).  Lake’s  [3]  dimples  were 
drilled  with  a  5.08  cm  ball  end  mill  to  a  depth  of  1 .588  mm.  Ultimately,  they  took 
elliptical  shape  due  to  the  curvature  of  the  blade  with  a  chord-wise  length  of  15.13  mm 
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and  a  span-wise  length  of  17.53  mm  with  a  center-to-center  spacing  of  2.22  cm.  Lake  [3] 
installed  dimples  at  three  axial  chord  locations  (50%,  55%,  and  65%)  and  researched 
them  using  a  linear  cascade  of  eight  blades  with  the  dimpled  blade  in  the  fifth  position. 

He  tested  operating  conditions  with  Reynolds  numbers  of  50,000,  100,000,  and  200,000. 
Each  case  led  to  a  loss  reduction  for  the  blade. 

For  his  research,  Rouser  [1]  used  the  same  blade  as  Lake  [3].  He  repeated  each  of 
Lake’s  [3]  cases  down  to  Reynolds  numbers  of  25,000.  He  also  experimented  with 
asymmetric  dimples  by  filling  half  of  each  dimple  with  Playdoh®  in  the  span-wise 
direction.  Ultimately,  he  concluded  that  the  full  dimples  at  65%  of  axial  chord  performed 
most  efficiently. 

Casey  [4]  also  used  the  same  blade  as  Lake  [3]  and  Rouser  [1].  In  order  to  test 
the  effectiveness  of  multiple  rows  of  dimples,  he  drilled  another  row  of  dimples  identical 
to  the  row  at  65%  of  axial  chord  at  76%  of  axial  chord.  He  also  looked  at  a  single  row  of 
dimples  at  65%  of  axial  chord  with  4.44  cm  spacing  rather  than  the  2.22  cm  spacing 
previously  used  by  Lake  [3]  and  Rouser  [1].  He  found  no  appreciable  differences 
between  the  different  dimple  patterns  and  concluded  that,  for  ease  of  manufacturing,  the 
single  row  at  65%  of  axial  chord  with  a  4.44  cm  dimple  spacing  was  the  best  choice  for 
boundary  layer  control. 


Figure  9;  Experimental  blade  after  modification  by  Casey  with  dimples  at  50%,  55%, 
65%,  and  76%  of  axial  chord 
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CHAPTER  3:  COMPUTATIONAL  STUDY 


Originally,  Rouser  [1]  created  a  eomputational  model  for  an  unmodified  Pak-B 
blade.  To  do  so,  he  eonstrueted  a  two-dimensional  grid  from  a  set  of  eoordinates  known 
to  define  the  Pak-B  profile.  Then,  the  two-dimensional  grid  was  extruded  span-wise  to 
ereate  a  three-dimensional  grid  model  of  the  basie  Pak-B  shape.  His  grid  ineluded 
struetured  grids  on  the  blade  surfaee  and  within  the  boundary  layer  and  unstruetured  grids 
in  the  rest  of  the  blade  passage.  The  baseline  grid  developed  by  Rouser  [1]  was  modified 
to  ereate  the  grid  used  in  this  researeh.  The  grid  topology  is  diseussed  in  greater  detail  in 
seetion  3.2. 

The  grid  model  was  simulated  with  an  axial  inlet  Reynolds  number  of  25,000.  A 
time-aeeurate  solution  was  developed  using  the  detaehed  eddy  simulation  turbulenee 
model. 

The  purpose  of  this  researeh  is  to  run  several  eases  using  the  detaehed  eddy 
simulation  turbulenee  model  and  determine  the  ability  of  the  model  to  make 
eomputations  for  use  in  predieting  future  results  from  experiments  for  this  applieation. 
Also,  the  results  may  enhanee  knowledge  of  the  physies  surrounding  flows  around  low 
pressure  turbine  blades  and  within  dimples  on  these  blades.  Unfortunately,  the  three- 
dimensional  grids  used  for  this  researeh  do  not  represent  aetual  physieal  struetures.  They 
negleet  the  effeets  of  blade  eharaeteristies  sueh  as  taper,  twist,  and  the  blade’s  rotating 
motion.  Therefore,  the  CFD  analysis  will  aeeurately  model  the  flow  field  only  on  the 
mid-span  of  the  blade.  Furthermore,  to  enhanee  the  two-dimensional  effeets,  the  blade  is 
modeled  with  an  infinite  span  and  the  easeade  is  modeled  with  an  infinite  number  of 
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blades.  Ultimately,  the  primary  reason  for  this  researeh  is  to  better  visualize  the  flow 
field  and  flow  struetures  within  the  field. 

3.1  Computational  Software 

Three  eommereial  software  paekages  were  used  during  the  eourse  of  this  researeh; 
Gridgen®  vl5,  Fluent®  v6.1.22,  and  Fieldview®  v9.  Gridgen®  is  a  grid  generation 
software  paekage  designed  to  interfaee  with  various  CFD  paekages.  The  researeh  grids 
were  ereated  and  modified  in  Gridgen®.  This  ean  be  a  delieate  proeess  sinee  the 
software  paekages  will  not  aeeurately  model  the  flow  field  without  a  proper  grid.  Eaeh 
node  within  the  grid  represents  a  point  in  spaee  at  whieh  the  governing  equations  are 
solved.  Grids  are  ereated  by  properly  plaeing  these  nodes  within  the  various  regions  that 
make  up  the  grid.  If  the  grid  points  are  plaeed  too  far  apart,  the  solver  eannot  aeeurately 
resolve  the  flow  field.  On  the  other  hand,  if  the  grid  points  are  too  near  eaeh  other,  the 
number  of  eells  inereases  dramatieally  and  a  prohibitive  amount  of  memory  and 
eomputing  time  is  neeessary.  There  are  two  types  of  grids;  struetured  and  unstruetured. 
On  struetured  grids,  the  number  of  grid  points  on  opposing  boundaries  must  be  the  same 
and  a  grid  point  on  one  boundary  eorresponds  to  a  grid  point  on  another  boundary. 
Neither  of  these  requirements  must  be  met  for  unstruetured  grids. 

Gridgen®  was  also  used  to  set  boundary  eonditions  for  the  grid  models.  Eaeh 
ease  used  a  veloeity  inlet  boundary  eondition  where  the  flow  enters  the  blade  passage  and 
a  pressure  outlet  eondition  where  the  flow  exits.  Also,  both  the  piteh-wise  and  span-wise 
faees  of  the  grid  were  set  with  periodie  boundary  eonditions  in  order  to  simulate  the 
infinite  linear  easeade  and  the  infinite  length  of  the  blade,  respeetively.  The  span-wise 
periodie  boundary  eondition  simulates  an  infinite  row  of  dimples  on  the  blade. 
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After  the  grid  was  ereated,  it  was  saved  as  a  Fluent®  ease  file  from  whieh  the  flow 
field  ean  be  resolves.  Fluent®  uses  a  finite  volume  method  to  diseretize  and  solve  the 
Navier-Stokes  equations  for  the  flow  field.  Finite  volume  methods  are  derived  from  the 
integral  form  of  the  equations  of  motion  and  then  solved  by  the  applieation  of  the 
eonservations  laws.  Eaeh  eell  within  the  grid  makes  up  a  finite  volume  in  whieh  all  flow 
properties  are  assumed  to  be  eonstant.  In  eaeh  finite  volume,  physieal  quantities  sueh  as 
mass  must  be  eonserved.  Fluent®  makes  these  ealeulations  by  integrating  density  in  eaeh 
eell  to  determine  its  instantaneous  mass.  The  time  rate  of  ehange  of  this  ealeulated  mass 
equals  the  mass  flux  through  the  eell.  These  rates  are  then  mathematieally  analyzed 
based  on  the  properties  in  adjaeent  eells  to  determine  flow  properties  of  eaeh  of  the  eells 
in  the  domain. 

Solutions  were  eomputed  with  the  following  boundary  eonditions.  The  inlet 
veloeity  eomponents  were  1.77804  m/s  in  the  x-direetion  and  1.221704  m/s  in  the  y- 
direetion.  This  eorresponds  to  an  inlet  angle  of  34.5'’.  The  inlet  pressure  was  set  to 
101,325  Pa  and  the  outlet  pressure  to  -4.4  Pa.  The  no-slip  eondition  was  set  at  the  blade 
surfaee.  The  DBS  model  eonstant  was  set  at  the  default  value,  0.65. 

Eaeh  grid  utilized  Eluent’s®  unsteady  solvers  to  ereate  a  time-aoeurate  solution. 

Eor  the  unsteady  solver,  the  time  step  is  set  to  a  eonstant  value  for  every  eell  within  the 
grid.  Though  this  leads  to  a  more  aeeurate  transient  solution  than  the  steady  solvers  ean 
provide,  it  also  requires  that  the  time  step  be  set  to  a  level  low  enough  to  resolve  the  most 
time-sensitive  region  of  the  flow  field.  The  unsteady  solver  eomputes  as  many  time  steps 
as  the  user  requires  with  a  set  maximum  number  of  iterations  per  time  step.  If  the 
eonvergenee  eriteria  set  by  the  user  are  met  during  the  eomputations  for  any  partieular 
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time  step,  the  partieular  time  step  cuts  short  and  the  solver  proceeds  to  the  next  time  step. 
For  this  research,  the  time  step  was  set  to  .00005  s  and  the  convergence  criteria  was  10'^. 
The  cases  were  run  on  the  Compaq®  CS40/45  at  the  Mission  Shared  Resources  Center 
on  between  24-32  processors  per  case. 

Post  processing  operations  were  performed  in  Fieldview®,  which  is  more 
functional  than  Fluent®  in  that  regard.  While  running  the  solutions  in  Fluent®,  the  files 
were  set  up  to  write  out  Fieldview®  case  and  data  fdes  every  .00075  s  of  flow  time, 
which  corresponds  to  15  time  steps.  This  allows  the  user  to  analyze  the  flow  field  at  a 
particular  instant  in  time  or  over  many  time  steps  during  the  run.  Fieldview’s® 
capabilities  include  calculations  of  the  separation  and  reattachment  locations,  velocity 
contours,  and  visualization  of  the  flow  field  using  streamlines. 

Fieldview®  was  also  used  to  make  calculations  and  show  diagrams  of  the  boundary 
layer.  To  do  so,  the  velocity  magnitude  was  plotted  along  a  2-cm  vector  normal  to  the 
surface  of  the  blade  at  its  mid-span.  The  boundary  layer  measurements  were  taken  at 
67.2%,  73.0%,  79.3%,  84.8%,  and  89.8%  of  the  axial  chord. 

3.2  Grid  Characteristics 

Initially,  a  file  containing  the  coordinates  that  describe  a  Pak-B  profile  was 
imported  into  Gridgen®.  This  profile  was  extruded  in  a  span-wise  direction  to  create  a 
three-dimensional  blade  with  the  Pak-B  shape  and  a  span  of  1 1 .28%  of  axial  chord.  Both 
grids  were  composed  of  both  structured  and  unstructured  regions.  In  both  cases,  the 
surfaces  of  the  blade  and  the  boundary  layer  around  the  blade  were  constructed  using  a 
structured  grid.  These  models  included  quadrilateral  cells  on  the  surface  transitioning 
into  hexahedral  cells  within  the  boundary  layer  region.  As  the  distance  from  the  blade 
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surface  increased,  the  grid  model  transitioned  to  an  unstructured  format  filled  with 
tetrahedral  cells.  For  the  dimpled  grid  model,  the  dimple  surface  was  modeled  by 
triangular  cells.  When  these  were  extruded  into  the  unstructured  region  above  the 
dimple,  the  cells  became  prisms. 

3,2,1  Grid  Topology  for  Unmodified  Pak-B  Blade 

For  the  unmodified  blade,  the  grid  that  Rouser  [1]  created  in  his  research  was 
modified  for  the  DBS  turbulence  model.  The  solution  was  resolved  on  the  Compaq® 
CS40/45  at  the  MSRC  on  24-32  processors  to  a  flow  time  of  .35  s.  On  this  grid,  5 1  nodes 
were  placed  across  the  span.  To  better  resolve  flow  structures,  a  higher  cell  density  was 
desired  in  regions  where  high  flow  gradients  were  expected.  Therefore,  more  cells  were 
located  near  the  leading  edge  and  the  trailing  edge  and  near  the  separation  location  and 
over  the  suction  surface  of  the  blade.  The  Pak-B  shape  is  defined  by  517  nodes.  The 
grids  on  the  blade  surface  were  then  extruded  normal  to  the  blade  to  create  the  boundary 
layer  region.  This  region  includes  5 1  nodes  on  the  normal  vector  starting  with  a  spacing 
of  .000080356  m.  The  spacing  grows  as  the  distance  from  the  blade  increases  with  a 
growth  rate  of  1 . 1 .  The  surface  of  this  domain  is  then  extruded  again  to  enclose  the 
remainder  of  the  expected  turbulent  region  as  required  by  the  DBS  model.  This  region 
includes  16  nodes  along  the  vector  normal  to  the  blade  and  has  an  initial  spacing  of 
.010178  m.  The  spacing  increases  with  a  growth  rate  of  1.25.  The  outer  boundary  of  the 
outlying  unstructured  region  includes  15  nodes  pitch- wise,  5  nodes  span- wise,  and  231 
nodes  stream- wise.  Ultimately,  the  grid  for  the  unmodified  blade  includes  1,801,522 
nodes  and  3,897,601  cells.  These  values  represent  an  increase  of  about  1  million  nodes 
and  2.5  million  cells  over  the  grid  originally  created  by  Rouser  [1].  The  primary 
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modifications  include  doubling  the  thickness  of  the  struetured  region  near  the  wall  to 
ensure  that  it  eaptured  the  turbulent  structures  in  the  flow  field  and  inereasing  the  number 
of  grid  points  in  the  boundary  layer  region  near  the  wall  to  more  accurately  resolve  flow 
field  eharaeteristies  in  this  region. 

3,2,2  Grid  Topology  for  Pak-B  Blade  with  Dimple  at  65%  of  Axial  Chord 

For  the  dimpled  blade,  the  grid  that  Casey  [4]  ereated  in  his  research  was  modified 
for  the  DBS  turbulence  model.  The  solution  was  resolved  on  the  Compaq®  CS40/45  at 
the  MSRC  on  24-32  processors.  Casey  [4]  created  a  spherical  database  structure  with  a 
diameter  of  5.08  cm  and  merged  it  with  the  database  for  the  unmodified  blade.  By 
cutting  the  sphere  at  the  intersection  between  its  surface  and  the  blade  surface  and 
ereating  a  domain  on  the  surface  of  the  spherical  database  structure,  Casey  [4]  ereated  a 
domain  on  the  dimple  centered  at  65%  of  axial  chord  and  with  a  depth  of  1 .588  mm.  The 
line  of  interseetion  between  the  blade  surfaee  and  the  sphere  eontained  360  equally 
spaced  nodes.  Above  the  dimple,  an  unstructured  region  was  included  to  the  same  height 
as  the  boundary  layer  region  around  the  rest  of  the  blade.  The  dimple  shape  was  extruded 
through  this  region  so  that  the  top  face  included  it.  Above  this  region,  another 
unstruetured  bloek  was  created  to  the  same  height  that  the  turbulent  region  was  modeled 
above  the  rest  of  the  block.  This  region  did  not  include  the  dimpled  shape.  Ultimately, 
the  grid  for  the  dimpled  blade  included  2,899,416  nodes  and  5,020,832  cells.  These 
values  represent  an  increase  of  about  1  million  nodes  and  2  million  cells  over  Casey’s  [4] 
grid.  The  primary  modifications  include  doubling  the  thiekness  of  both  the  struetured 
region  near  the  wall  and  the  unstructured  region  above  the  dimple  in  order  to  ensure  that 
it  eaptured  the  turbulent  structures  of  the  flow  field. 
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Fig.  10  illustrates  the  basie  grid  for  the  mid-span  of  the  blade  passage  for  eaeh  ease. 


Fig.  1 1  illustrates  the  basie  grid  in  the  boundary  layer  for  eaeh  ease. 
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Figure  10:  Grid  topology  ofPak-B  mid-span  passage  [Rouser,!] 


Figure  11:  Grid  topology  of  boundary  layer 
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3.3  Computational  Results  for  Pak-B  Unmodified  Blade 
3.3,1  Separation  Location 

Fig.  12  shows  a  diagram  of  the  locations  of  separation  and  reattachment  as 
calculated  by  Fieldview®.  To  provide  a  sense  of  perspective,  the  green  dots  are  located 
at  67.2%,  73.0%,  79.3%,  84.8%,  and  89.8%  of  axial  chord.  The  diagram  indicates  that 
no  separation  occurs  until  the  flow  field  has  passed  89.8%  of  axial  chord.  It  also 
indicates  that  the  flow  attempts  to  reattach  to  the  blade  surface  at  many  positions 
although  there  are  no  reattachment  locations  that  totally  the  blade  surface.  Rather,  the 
reattachment  lines  follow  the  direction  of  the  flow.  For  comparison,  Casey  [4]  found  that 
Fieldview®  calculated  the  location  of  separation  for  an  unmodified  blade  at  between 
64.5%  and  68.5%. 


■  separation 

■  reattachment 


Flow  Direction  - ^ 


Figure  12;  Separation  and  reattachment  locations  for  unmodified  Pak-B  blade  run  at  Re 
25k  using  DBS  turbulence  model  at  0.3495  s 
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3,3,2  Velocity  Profile 


Fig.  13  shows  the  velocity  magnitude  contours  of  the  computational  flow  field  for 
an  unmodified  Pak-B  blade  on  a  plane  located  at  the  mid-span  of  the  blade.  It  illustrates 
the  presence  of  vortical  structures  within  the  turbulent  region  of  the  flow  field  as 
indicated  by  the  increase  in  velocity  contained  within  the  separated  region.  Also,  Fig.  1 1 
shows  that  the  flow  separates  from  the  blade  surface  at  the  position  where  the  boundary 
layer  begins  to  move  away  from  the  blade.  Note  that  this  position  is  much  farther 
upstream  than  the  position  that  Fieldview®  calculated  for  separation.  The  separated 
region  of  the  flow  field  steadily  increases  in  size  beginning  at  the  location  of  primary 
separation.  Though  the  rate  of  growth  slows  as  the  flow  moves  downstream,  the 
separated  region  appears  to  grow  in  size  until  the  flow  reaches  the  trailing  edge  of  the 
blade.  Span-wise  uniformity  of  the  contour  plots  were  checked  and  no  significant 
differences  were  discovered. 

P  1.8^6 


“  0.000 


Figure  13:  Contour  plot  of  normalized  velocity  magnitude  for  unmodified  Pak-B  blade  at 
Re  25k  using  DBS  turbulence  model  at  0.3495  s 
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Fig.  14  shows  a  close-up  of  the  veetor  field  that  defines  the  turbulent  region  of  the 
flow  on  a  plane  loeated  at  the  mid-span  of  the  blade.  A  large  vortieal  strueture  shedding 
from  the  blade  dominates  this  region  of  the  flow  field.  Fig.  15  depiets  another  elose-up 
of  the  veetor  field  nearer  the  blade  surfaee.  The  red  lines  loeated  beneath  the  veetor  field 


illustrate  that  the  flow  attempts  to  reattaeh  to  the  blade  surface  within  the  turbulent 
region.  The  large  baekflow  present  within  the  turbulent  region  of  the  flow  field  rejeets 
the  attempts  of  the  boundary  layer  to  reattaeh.  It  is  likely  that  the  previously  noted 
vortieal  structure  has  a  large  effeet  on  these  reattachment  attempts  sinee,  as  evideneed  by 
the  greater  density  of  reattaehment  lines,  the  reattachment  oeeurs  more  strongly 
immediately  after  the  vortieal  strueture  passes  downstream.  These  results  match  previous 


researeh  performed  by  Flalstead  et  al.  [3]  and  Casey  [4]. 


.  \ 


V  vAa  \  \ 


I 


'  .  V.  \.  V.  \,\  A;A  .A,a  .\  \  A 
t  A  A  . A  A  A  A 

W\  A  .A  A.  \ ,  A  \  A  A 

■  A  V  A  ,  A  A  . '  A  A  A  V 

'a'  a  a'  a'  a  a  a  a  a  V  A  , 

A  A  A  A  A  '  '  ' 


•  ■  ■  - 

•  '  ^  W  A-. 


A  A  i  V  ^  -A 


Figure  14:  Veetor  plot  of  velocity  magnitude  in  turbulent  region  for  unmodified  Pak-B 
blade  at  Re  25k  using  DBS  turbulenee  model  at  0.3495  s 
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Figure  15;  Vector  plot  of  velocity  magnitude  in  turbulent  region  for  unmodified  Pak-B 
blade  at  25k  using  DBS  turbulence  model  at  0.3495  s 

3,3,3  Streamline  Visualizations 

Fig.  16  shows  a  side  view  of  a  set  of  surface  streamlines  along  the  entire  span  of 
the  blade.  The  streamlines  show  definite  separation  from  the  blade  surface  between  65% 
and  70%  of  axial  chord. 
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Figure  16;  Streamline  visualization  of  flow  field  for  unmodified  Pak-B  blade  at  Re  25k 
using  DBS  turbulence  model  at  0.3495  s  of  flow  time 

3.3,4  Computational  Boundary  Layer  Measurements 

Fig.  17  shows  profiles  of  the  boundary  layer  constructed  for  the  unmodified  Pak-B 
computational  model  along  the  mid-span  of  the  blade.  At  67.2%  of  axial  chord,  there  is  a 
very  slight  inflection  point  located  near  the  blade  surface.  That  inflection  point  indicates 
that  the  flow  has  already  begun  to  separate  from  the  blade  surface  at  this  location.  At 
73.0  %  of  axial  chord,  there  is  a  small  separation  bubble  forming  in  the  flow  field.  As  the 
flow  moves  downstream,  the  separation  bubble  steadily  grows  in  size. 
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Figure  17;  Boundary  layer  profiles  for  unmodified  Pak-B  blade  at  Re  25k  using  DBS 
turbulence  model  at  0.3495  s 


For  purposes  of  comparison,  Fig.  18  shows  boundary  layer  profiles  determined 
both  computationally  and  experimentally  by  Casey  [4]  for  an  unmodified  blade  using  a 
laminar  solver. 
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Figure  18;  Boundary  layer  comparisons  between  laminar  CFD  model  and  experiment  for 
an  unmodified  Pak-B  blade  at  Re  25k  and  0%  turbulence  [Casey ,4] 

3.4  Summary  and  Conclusions  for  Unmodified  Blade 

According  to  calculations  made  by  Fieldview®,  the  flow  does  not  separate  to  any 
significant  degree  until  it  has  reached  approximately  90%  of  axial  chord.  This  result  does 
not  match  either  the  velocity  magnitude  contour  plot  in  Fig.  13  or  the  boundary  layer 
profiles  determined  for  this  research  shown  in  Fig.  17  or  determined  by  Casey  [4]  in  his 
research  shown  in  Fig.  18.  Each  of  these  other  data  sources  indicates  that  separation 
begins  to  occur  closer  to  65%  of  axial  chord.  The  difference  between  this  data  and  these 
other  data  sources  indicates  that  the  separation  location  calculated  by  Fieldview®  is 
incorrect  and  should  be  disregarded.  As  such,  it  was  not  possible  to  determine  an  exact 
separation  location.  Flowever,  analysis  of  the  other  data  sources  indicates  that  separation 
occurred  near  the  position  discovered  by  Casey  [4]. 

The  velocity  magnitude  contour  plot  and  the  associated  vector  plots  provide  much 
more  reliable  interpretation  of  the  data.  According  to  the  contour  plot,  the  flow  separates 
from  the  blade  surface  around  65%  of  axial  chord  and  does  not  reattach.  Also,  it  shows 
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that  the  flow  fleld  in  the  turbulent  region  is  dominated  by  vortieal  struetures  that  induee  a 
signifieant  baek  flow.  These  results  elosely  mateh  similar  results  reported  by  Casey  [4], 
These  images  verify  that  the  DBS  model  ean  resolve  turbulent  struetures  in  the  flow  field. 
This  makes  sense  sinee  the  DBS  model  is  equivalent  to  the  BBS  model  in  any  region  with 
signifieant  amounts  of  turbulenee. 

The  boundary  layer  profiles  are  primarily  used  to  determine  the  thiekness  of  the 
boundary  layer.  In  this  ease,  the  boundary  layer  signifieantly  inereases  in  thiekness  as  the 
flow  moves  downstream.  Though  these  results  do  not  perfeetly  mateh  Casey’s  [4] 
experimental  results,  they  do  eompare  favorably  to  his  eomputational  results,  for  whieh 
he  employed  a  laminar  solver.  At  67.2%  of  axial  ehord,  both  sets  of  eomputational  data 
mateh  eaeh  other  elosely.  However,  as  the  flow  moves  downstream,  the  results  from  this 
researeh  show  a  smaller  boundary  layer  than  the  results  from  Casey’s  [4]  eomputational 
researeh.  This  lower  reeorded  thickness  for  the  DBS  model  more  closely  matches  the 
experimentally  observed  boundary  layer  thickness  recorded  by  Casey  [4].  This 
represents  the  first  significant  advantage  of  the  DBS  turbulence  model  over  the  laminar 
solver.  It  should  be  noted  that  Casey’s  [4]  results  are  time-averaged  whereas  the  results 
from  this  research  are  not.  However,  this  difference  is  not  expected  to  have  a  great  effect 
on  this  comparison. 

Bor  most  of  the  research  data,  the  DBS  turbulence  model  proved  either  as  accurate 
or  more  accurate  than  the  laminar  solver  for  an  unmodified  Pak-B  blade.  Since  Casey  [4] 
already  determined  that  the  laminar  solver  is  more  accurate  than  other  currently  available 
turbulence  models,  future  researchers  should  strongly  consider  using  the  DBS  model  in 
their  research. 
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3.4  Results  for  Pak-B  Blade  with  Dimple  at  65%  of  Axial  Chord 

3,4,1  Separation  Location 

Fig.  19  shows  a  diagram  of  the  locations  of  separation  and  reattachment  as 
calculated  by  Fieldview®  at  0.15  s  of  flow  time.  The  diagram  shows  that  the  flow 
separated  from  the  blade  as  it  enters  the  dimple  and  reattaches  as  it  exits  the  dimple. 
Subsequently,  the  flow  separates  again  at  approximately  72%  of  axial  chord.  Beyond  this, 
there  is  no  significant  line  of  reattachment.  For  comparison,  Casey  [4]  found  the 
contained  vortex  within  the  dimple  and  calculated  the  separation  line  at  71%  of  axial 
chord.  Subsequently,  he  discovered  that  the  flow  reattaches  at  83%  of  axial  chord.  Even 
though  the  results  from  this  study  show  separation  downstream  of  Casey’s  results,  Casey 
also  discovered  that  the  separation  location  moves  upstream  as  the  flow  time  progressed 
so  the  results  are  consistent. 


■  Separation 

■  Reattachment 
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Figure  19;  Separation  and  reattachment  lines  for  Pak-B  blade  with  dimple  at  65%  of  axial 
chord  for  Re  25k  using  DES  turbulence  model  at  0.15  s 

3,4,2  Velocity  Field  Plots 

Eig.  20  shows  the  velocity  magnitude  contours  of  the  computational  flow  field  for 
a  Pak-B  blade  with  a  dimple  located  at  65%  of  axial  chord  on  a  plane  located  at  the  mid¬ 
span  of  the  blade.  It  illustrates  the  presence  of  vortical  structures  within  the  turbulent 
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region  of  the  flow  field.  Also,  Fig.  20  shows  that  the  flow  separates  from  the  blade 
surface  at  the  position  where  the  boundary  layer  begins  to  move  away  from  the  blade. 
Note  that  the  presence  of  the  dimple  causes  separation  to  occur  farther  downstream  than 
for  the  unmodified  case.  Span-wise  uniformity  of  the  contour  plots  were  checked  and  no 
significant  differences  were  discovered. 


Figure  20:  Normalized  velocity  magnitude  contour  plot  for  Pak-B  blade  with  dimple  at 
65%  of  axial  chord  at  Re  25k  using  DBS  turbulence  model  at  0.15  s 

Fig.  21  shows  the  vector  field  for  the  turbulent  region  of  the  dimpled  case.  Note 
the  presence  of  vortical  structures  though  the  structure  is  much  closer  to  the  trailing  edge 
of  the  blade  than  for  the  unmodified  case. 
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Figure  10;  Vector  field  of  velocity  magnitude  in  turbulent  region  for  Pak-B  blade  with 
dimple  at  65%  of  axial  chord  at  Re  25k  using  DBS  turbulence  model  at  0.15  s 


Fig.  22  shows  the  vector  field  of  the  flow  as  it  enters  the  dimple.  Note  that  the 
flow  travels  uniformly  downstream  until  it  enters  the  dimple  at  which  point  it  gets  caught 
up  in  the  dimple’s  contained  vortex. 
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Figure  22;  Vector  field  of  velocity  magnitude  of  flow  entering  dimple  for  Pak-B  blade 
with  dimple  at  65%  of  axial  chord  at  Re  25k  using  DBS  turbulence  model  at  0.15  s 

3,4,3  Streamline  Visualizations 

Fig.  23  displays  the  behavior  of  the  streamlines  of  the  flow  field  released  from  a 
row  of  seeds  located  at  50%  of  axial  chord.  The  streamline  visualization  shows  that  the 
vortex  contained  within  the  dimple  captures  the  streamlines  and  pulls  them  into  the 
dimple.  Subsequently,  the  dimple  expels  the  flow  down  the  centerline.  It  stays  at  the 
centerline  until  it  enters  the  turbulent  region  of  the  flow  field  and  the  flow  becomes  very 
chaotic. 
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Figure  23;  Streamline  visualization  of  flow  field  for  Pak-B  blade  with  dimple  at  65%  of 
axial  ehord  at  Re  25k  using  DBS  turbulence  model  at  0.15  s 

Fig.  24  shows  a  side  view  of  the  streamlines  within  the  dimple.  The  streamlines 
clearly  show  the  contained  vortex  within  the  dimple. 


Figure  24:  Streamline  visualization  of  dimple  within  flow  field  of  Pak-B  blade  with 
dimple  at  65%  of  axial  chord  at  Re  25k  using  DBS  turbulence  model  at  0.15  s 
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3,4,4  Boundary  Layer  Profiles 


Fig.  25  shows  profiles  of  the  boundary  layer  constructed  for  the  Pak-B 
computational  model  with  a  dimple  located  at  65%  of  axial  chord  along  the  mid-span  of 
the  blade.  At  67.2%  and  73.0%  of  axial  chord,  there  is  little  evidence  of  separation.  By 
79.3%  of  axial  chord,  there  is  a  small  separation  bubble.  As  the  flow  moves  downstream, 
the  separation  bubble  steadily  grows  in  size. 


Normalized  Velocity  Magnitude  (Umeasured/Uiniet) 

Figure  25:  Boundary  layer  profiles  for  Pak-B  blade  with  dimple  at  65%  of  axial  chord  at 
Re  25k  using  DBS  turbulence  model  at  0.15  s 

For  purposes  of  comparison.  Fig.  26  shows  Casey’s  experimental  and 
computational  results  for  this  case. 
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Figure  26;  Boundary  layer  eomparisons  between  CFD  and  experiment  for  a  Pak-B  blade 
with  dimple  at  65%  of  axial  ehord  at  Re  25k  and  1%  Tu  [Casey, 4] 

3.5  Summary  and  Conclusions  for  Pak-B  Blade  with  Dimple  at  65%  of 
Axial  Chord 

According  to  calculations  made  by  Fieldview®,  the  flow  separates  around  73%  of 
axial  chord.  This  result  matches  other  data  sources  such  as  the  velocity  magnitude  plots 
and  boundary  layer  profiles.  This  differs  from  the  results  for  the  unmodified  blade.  This 
difference  may  result  from  the  major  differences  in  the  respective  flow  fields.  For  this 
case,  the  dimple  works  to  change  the  flow  field  from  laminar  to  turbulent.  Flowever,  with 
the  unmodified  blade,  there  is  no  structure  designed  to  accomplish  this  so  the  solver  must 
decide  where  the  flow  transitions.  On  the  other  hand,  these  results  do  not  show  the 
subsequent  reattachment  discovered  both  experimentally  and  computationally  by  Casey 
[4].  However,  since  the  flow  field  did  not  develop  as  greatly  for  this  research,  it  is 
possible  that  further  development  would  show  reattachment. 

According  to  the  contour  plot,  the  flow  separates  from  the  blade  surface  well 
downstream  of  the  dimple.  Also,  it  shows  that  the  flow  field  in  the  turbulent  region  is 
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marked  by  vortical  structures  that  induce  a  significant  back  flow.  However,  these  results 
do  not  show  as  much  turbulent  activity  in  the  form  of  localized  regions  of  high  velocity 
as  reported  by  Casey  [4]  in  his  computational  results.  Again,  this  difference  may  result 
from  the  fact  that  Casey  [4]  was  able  to  develop  the  flow  field  to  a  greater  extent  than 
reached  in  this  research.  These  images  also  verify  that  the  DBS  model  can  resolve 
turbulent  structures  in  the  flow  field.  This  makes  sense  since  the  DBS  model  is 
equivalent  to  the  BBS  model  in  any  region  with  significant  amounts  of  turbulence. 

The  streamlines  show  the  same  behavior  as  reported  by  Casey  [4]  in  his  research. 
Specifically,  they  show  that  the  dimple  vortex  captures  the  streamlines  and  pulls  them 
into  the  dimple  and  expels  them  along  the  centerline  of  the  dimple. 

The  boundary  layer  profiles  are  primarily  used  to  determine  the  thickness  of  the 
boundary  layer.  In  this  case,  the  boundary  layer  increases  in  thickness  as  the  flow  moves 
downstream.  Though  these  results  do  not  perfectly  match  Casey’s  [4]  experimental 
results,  they  do  compare  favorably  to  his  computational  results,  for  which  he  employed  a 
laminar  solver.  At  67.2%  of  axial  chord,  both  sets  of  computational  data  match  each 
other  closely.  At  73.0%  of  axial  chord,  Casey’s  [4]  computational  data  shows  a  small 
separation  bubble  whereas  neither  his  computational  data  nor  the  data  from  this  research 
shows  this  feature.  As  the  flow  moves  downstream,  the  results  from  this  research  match 
Casey’s  [4]  experimental  research  better  than  the  results  from  Casey’s  [4]  computational 
research,  ft  should  be  noted  that  Casey’s  [4]  results  are  time-averaged  whereas  the 
results  from  this  research  are  not.  However,  this  difference  is  not  expected  to  have  a 
great  effect  on  this  comparison. 
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For  most  of  the  research  data,  the  DBS  turbulence  model  proved  either  as  accurate 
or  more  accurate  than  the  laminar  solver  for  a  Pak-B  blade  with  a  dimple  located  at  65% 
of  axial  chord.  These  results  strengthen  the  contention  that  future  researcher  should 
strongly  consider  the  DBS  turbulence  model. 
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Chapter  4:  Summary  and  Recommendations 


As  flight  technology  has  advanced,  the  United  States  military  has  begun  to  rely 
more  heavily  on  unmanned  aerial  vehicles,  such  as  the  Navy’s  X-47  Pegasus,  to  complete 
essential  missions.  These  aircraft  often  fly  at  high  altitudes  and  low  speeds,  conditions 
that  lead  to  low  Reynolds  numbers.  Ultimately,  these  conditions  cause  flow  separation 
off  the  blades  of  low  pressure  turbines  and  thereby  greatly  decrease  the  effectiveness  of 
the  aircraft  so  that  variables  such  as  range  and  endurance  are  decreased. 

Given  these  difficulties,  many  resources  have  been  invested  into  research  of  both 
active  and  passive  separation  control  methods.  Advantageously,  passive  controls  help  to 
conserve  the  aircraft’s  energy  and  are  particularly  useful  in  systems  in  which  flow 
conditions  maintain  fairly  steady  values.  Researchers  have  shown  that  dimples  are 
among  the  more  successful  passive  control  devices.  Specifically,  spherical  dimples 
placed  at  65%  of  axial  chord  have  proven  to  optimize  the  loss  reduction  for  LPT  blades. 

Unfortunately,  the  CFD  research  has  lagged  behind  the  experimental  research.  To 
date,  no  turbulence  model  exists  that  can  accurately  resolve  the  entire  flow  field. 
Previously,  a  laminar  viscous  model  had  been  most  successful  at  solving  this  problem. 
The  current  research  aimed  to  investigate  the  effectiveness  of  the  detached  eddy 
simulation  turbulence  model.  Two  cases  were  studied  using  this  model:  an  unmodified 
Pak-B  blade  and  a  Pak-B  blade  with  a  spherical  dimple  located  at  65%  of  axial  chord. 

The  DBS  model  was  used  to  determine  the  separation  location,  prepare  velocity 
magnitude  contour  plots  and  vector  plots  of  the  flow  field,  visualize  the  streamlines 
associated  with  the  dimple,  and  create  boundary  layer  profiles  from  which  the  thickness 
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of  the  boundary  layer  eould  be  determined.  Overall,  the  DBS  turbulenee  model 
outperformed  the  laminar  solver. 

Regarding  the  use  of  the  DBS  model  to  analyze  dimples  on  low  pressure  turbine 
blades,  more  researeh  remains  to  be  done.  To  ensure  that  the  CBD  solutions  are  as 
physieally  aeeurate  as  possible,  rigorous  grid-  and  time-independenee  studies  should  be 
performed.  Also,  a  study  should  be  performed  to  either  assure  that  the  default  values  for 
the  turbulenee  model  optimize  its  aeeuraey  or,  if  not,  determine  what  values  would  do  so. 
Also,  future  researeh  should  endeavor  to  extend  the  amount  of  flow  time  for  whieh  a 
solution  is  obtained. 

In  eonelusion,  this  researeh  provided  information  regarding  the  ability  of  the  DBS 
turbulenee  model  to  aeeurately  resolve  the  flow  field  for  both  an  unmodified  Pak-B  blade 
and  a  similar  blade  with  a  dimple.  In  many  respeets,  the  DBS  model  proved  to  be  as 
sueeessful  as  any  other  solver  for  these  eases.  This  will  benefit  researehers  sinee  they  ean 
now  more  reliably  prediet  experimental  results  using  CBD  simulations. 
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Appendix  A:  Sample  Computer  Files 


Example  of  Fluent®  infile  used  in  researeh;  ealled  from  bateh  submission  seript; 

fde 

rod 

<lile  name  of  file  that  Fluent®  should  oall> 

q 

q 

file 

read-maoros 

maoros.som 

q 

q 

grid 

reorder 

reorder-domain 

q 

q 

q 

define 

models 

unsteady-2nd-order 

q 

q 

solve 

set 

time-step 

0.00005 

q 

dual-time-iterate 

2000 

20 

q 

q 

file 

wed 

<file  name  that  Fluent®  should  write  to> 

q 

q 

/parallel/timer/print 

q 

q 

exit 

y 
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Example  of  batch  submission  script  used  in  research; 

#!/bin/csh  -x 
#BSUB  -J  <job  name> 

#BSUB  -n  <number  of  processors> 

#BSUB  -W  <maximum  clock  time  desired  in  hh:mm  format> 

#BSUB  -q  Default 
#BSUB  -0  <output  fde> 

#BSUB  -e  <error  lile> 

#BSUB  -P  <16  character  account  number> 

cd  $WORK_DIR 

/usr/bin/rcp  <case  fde  path>  . 

/usr/bin/rcp  <data  fde  path>  . 

/usr/bin/rcp  <inlile  path>  . 

fluent  3ddp  -g  -ptmpi  -t  <number  of  processors>  -i  <infde  name>  >&  new.  out 

/usr/bin/rcp  <output  case  file  name>  ${msas}:<case  file  path>  . 

/usr/bin/rcp  <output  data  file  name>  ${msas};<data  file  path>  . 

/usr/bin/rcp  <out  file  name>  ${msas}:<data  file  path>  . 


Example  of  macro  file  used  in  research  to  periodically  write  Eieldview®  files;  called 
from  Eluent®  infile: 

(cx-macro-define 

‘((fv  .  “\nfile\nexport\nfieldview\n<file  name%t>\npressure\ndyn\ntotal\nx-vel\ny- 
vel\nz-vel\nq\n”) 

)) 
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